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Abstract 

The Neumann system on the 2-dimensional sphere is used as a 
tool to convey some ideas on the bi-Hamiltonian point of view on 
separation of variables. It is shown that, from this standpoint, its 
separation coordinates and its integrals of motion can be found in a 
systematic way. 

1 Introduction 

Separation of variables for the Hamilton- Jacobi equation is a classical topic 
that is still very much investigated and has connections with several impor- 
tant research fields, such as stationary reductions of soliton equations (see 
|8| and references therein), algebraic completely integrable systems [I], |TT| , 
Riemannian geometry || [lj^ ^(| |29| , Backlund transformations and Baxter's 
Q-operator KR |2T|, and string theory ftM. 



The aim of this paper is to present, in the simple example of the Neu- 



mann system, the main ideas of a new approach [Ty, |22], [HI to separation 
of variables. This point of view is based on the geometry of bi-Hamiltonian 
manifolds, and has been successfully applied to the stationary reductions of 
the KdV hierarchy [JD]] and to Toda lattices ||. 

The Neumann system is a well-known and very much studied mechanical 
system (see, e.g., [El, F2J, [15], 56]), given by a (mass 1) particle moving on 



the (unit) sphere under the influence of a quadratic potential V(x, y, z) 
\{a\x 2 + a 2 y 2 + a^z 3 ), where a\ < a 2 < (13 and (x,y,z) are Cartesian co- 
ordinates whose origin coincides with the center of the sphere. In 1859, 
Carl Neumann showed that the Hamilton- Jacobi equation of this system 
is additively separable in the so-called spheroconical (or elliptical spherical) 
coordinates (Ai,A 2 ), given by 

x 2 | y 2 | _ (A - Ai)(A - A 2 ) 

A — di A — CL2 A — (A — a x )(A — a 2 )(A — 0,3) ' 
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with di < Ai < a 2 < A 2 < a 3 . Once they have been introduced, one can 
check that they are separation coordinates, but the problem we want to 
address is to find out these coordinates in a systematic way. It is known that 
the separation variables of the Neumann system can be supplied by a Lax 
representation, but then the problem is to find such a representation for a 
given system, which is in general a quite difficult task. We will show that 
a careful study of the Neumann system will allow us to deduce in a natural 
way the spheroconical coordinates as separation variables. 

The plan of the paper is the following. In Section 2 we will introduce 
the idea (due to A. Nijenhuis) of the geometrization of a coordinate system 
by means of a suitable tensor field L. This allows one to find a kind of 
compatibility condition between L and a given Hamiltonian H, ensuring 
that the coordinates induced by L separate the Hamilton-Jacobi equation 
associated with H. This condition is an intrinsic (i.e., coordinate-free) form of 
the classical Levi-Civita separability conditions. Before applying this results 
to the Neumann system in Section 4, we will comment on their meaning in 
the theory of bi-Hamiltonian manifolds in Section 3. Section 5 is devoted to 
two additional results on the Neumann system, that can be easily obtained 
from our standpoint (using some facts highlighted in [I5||): A simple rule 
for the construction of the integrals of motion and the extension of the bi- 
Hamiltonian structure from the cotangent bundle T*S 2 of the sphere to a 
5-dimensional manifold, in order to give the Neumann system a bi-Hamil- 
tonian formulation, which is missing in T*S 2 . The final section contains some 
concluding remarks. 
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under the research project Geometry of Integrable Systems. 



2 The search for separation coordinates 

Let Q be an rz-dimensional manifold, that can be thought of as the configu- 
ration space of a mechanical system (the sphere, for the Neumann system), 
and let H be a function on T*Q (the Hamiltonian). In this section we will 
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present a method to look for the variables in which the Hamilton-Jacobi 
equation for H separates. This strategy will prove to be quite efficient in the 
case of the Neumann system. 

The first important point, due to Nijenhuis p5| , is the idea of the ge- 



ometrization of a coordinate system. If {q 1 } are local coordinates on an open 
subset U C Q, then the tensor field L of type (1,1) defined by 

L r = q r 

has vanishing Nijenhuis torsion, i.e., 

[LX, LY] - L[LX, Y] - L[X, LY] + L 2 [X, Y] = 

for every pair (X, Y) of vector fields on U. Viceversa, a tensor field L of type 
(1,1) whose torsion is zero induces local coordinates in a neighborhood of 
any point where the eigenvalues (Ai, . . . , A n ) of L are distinct. In many cases 
these eigenvalues are functionally independent, so that 

L JL = Xt ± 5 



i.e., they can be chosen as the coordinates associated with L. In the following 
we will suppose to be in this situation, that is, we will look for a tensor field 
L with vanishing torsion and functionally independent eigenvalues. 

Once we have replaced the coordinate system with the geometric object 
L, it is quite natural to look for a "compatibility condition" between L and 
a function H G C°°(T*Q) entailing that the coordinates given by L are 
separation variables for H. This condition can be obtained by lifting L from 
Q to T* Q by means of a procedure called "complete lifting" [[30], [(J . It gives 
rise to a torsionless tensor field iV of type (1, 1) on T*Q, described in fibered 
coordinates (q l ,Pi) as 

AT A = L i ± + ( d A _ H) Pk — , iV— = V— . 
dq l % dqi I dq l dq j J dpj ' dpi j dpj 

Now that we have H and iV on the same manifold T*Q, we can obtain the 
separability condition we are looking for as follows. First, we use H and iV to 
construct the 2-form ujh = d (N*dH), where iV* is the adjoint of N. Then we 
consider the Hamiltonian vector field Xh and the vector fields NXh, N 2 Xh, 
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. . . obtained by iteration, and we suppose that the distribution D H spanned 
by Xh, NXh, • • • , N n ~ x XH is n-dimensional (i.e., that these vector fields 
are linearly independent at any point). 

Theorem 1 In the above-mentioned hypotheses, the coordinates associated 
with L are separation variables for H if and only if the 2-form ujh annihilates 
the distribution Dm- 

u h \d h = . (1) 

The proof consists in writing the conditions 

u H (N k X H ,N l X H ) =0 , for k,l = 0,...,n- 1, (2) 

in the canonical coordinates (Ai, . . . , A„, //i, . . . , (J, n ), where the A« are the 
eigenvalues of L and the are their conjugate momenta. Then, taking into 
account that N is diagonal in these coordinates, it is not difficult to realize 
that the equations (fj) are equivalent to the classical Levi-Civita separability 
conditions ([|], p. 208). 

Therefore, we have found the separability condition between L and H. 
It is simply the vanishing of a suitable 2-form on the distribution generated 
from Xh by means of the complete lift N of L. We stress that this condition 
is a concise and, above all, intrinsic form of the Levi-Civita equations. This 
means that one can check the separability of a given Hamiltonian H in the 
coordinates associated with a tensor field L before computing these coordi- 
nates. Moreover, one can impose condition ([!]) on L to search for separation 
coordinates for H. In Section 4 we will show how to exploit this fact in order 
to systematically deduce the separability of the Neumann system. 



3 The bi- Hamiltonian meaning of the sepa- 
rability condition 

Before applying the results of the previous section to the Neumann system, 
let us make some comments on their geometrical meaning. 

1. Let Q, L, and iV be as in Section 2. Then the cotangent bundle of Q 
is a bi-Hamiltonian manifold [IS|. The first Poisson bracket is the canonical 
one, 

{F,G} = u(X F ,X G ) , 
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while the second one is given by 



{F,G}' = lo(NX f ,X g ) . 

If P (resp. P') is the Poisson tensor of {•,•} (resp. {•,•}'), this means that 
P' = NP. 

2. Assume that H fulfills the separability conditions (H). Then the distri- 
bution Dh is integrable, so that there (locally) exist n independent functions 
Hi, ... , H n which are constant on the leaves of Dh- They are (local) first in- 
tegrals for Xh- Notice however that the Hamiltonian vector field associated 
with the separable Hamiltonian H is not (in general, and in the particular 
example of the Neumann system) bi-Hamiltonian. 

3. The distribution Dh is bi-Lagrangian, that is, 

{H i} Hj} = , {Hi, Hj}' = for all i, j. 



As explained in [10 , |13|j , this condition characterizes the (integrable) Hamilto- 
nian systems that are separable in the coordinates induced by L. We remark 
in passing that bi-Lagrangian foliations play an important role in the study 



of special Kahler manifolds [16 



4. As we have already said in the previous section, the separation coor- 
dinates (Aj,/ij) are canonical for co (this is obvious) and for N, 

d d d d 

They are often called Darboux-Nijenhuis coordinates (since is sometimes 
called the Nijenhuis tensor) and have been used as separation coordinates in, 
e. g .,@||,0,i 



4 The case of the Neumann system 

In this section we will exploit the separability condition (|l|) to find the vari- 
ables of separation for the Neumann system. To this aim, we will look for a 
torsionless tensor field L of type (1, 1) on S 2 such that its complete lift iV on 
T*S 2 satisfies 

uh(X h ,NXh) = , (3) 
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where ujh = d(N*dH) and H is the Neumann Hamiltonian. Due to the 
form of the constraint and of the potential, it is quite natural to perform 
the computations in the local coordinates X = x 2 , Y = y 2 , parametrizing 
every connected component of the subset obtained by removing from S 2 its 
intersections with the coordinate planes. In these coodinates the Hamiltonian 
of the Neumann system has the form: 

H = 2 [X(l - X)p 2 x - 2XY PxPY + Y(l - Y)p 2 Y ] 

11 (4) 

+ o( a l - a 3) X + ^( a 2 - a 3) Y > 

where {px,Py) are the conjugate momenta of (X,Y). The unknown tensor 
field L can be written as 

L*(dX) = AdX + BdY 

L*(dY) = CdX + DdY , ^ ' 

where (A, B,C, D) are functions of the coordinates (X, Y) that must satisfy 
two additional conditions. The first one is that the torsion of L has to vanish, 
that is, 

B (C y - D x ) - BA X - DA Y + AA Y + CB Y = 
C (B x - A Y ) - AD X - CD Y + BC X + DD X = . 

The second one is the independency of the eigenvalues of L. Our aim is to 
find (A, B, C, D) in such a way to verify also the separability conditions (BD. 
Thus we need to compute the complete lift N of L, which turns out to be 
given by 



N*dX = AdX + BdY 
N*dY = CdX + DdY 



i\ ui = uoa -|- uai 

N*d Px = - \p x (B x - A Y ) + p Y (D x - C Y )] dY + Adp x + Cdp Y 
N*dp Y = [p x (B x - A Y ) + p Y (D x - C Y )} dX + Bdp x + Ddp Y . 



(6) 



Then (Q) becomes a differential equation where the coordinates (X,Y), the 
momenta (px,Py), and the unknown functions (A, B, C, D) and their deriva- 
tives appear polynomially. This suggests to seek for a solution which also 
depends polynomially on the coordinates. Trying the simplest solution, one 
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finds 



A = 


(03 


— ai)X + a\ 


B = 


(a 3 


-a 2 )X 


C = 


(03 


-oi)y 


D = 


(03 


- a 2 )F + a 2 



(7) 



(8) 



which leads to the tensor field L defined by 

L*(dX) = ai dX + Xd[(a 3 - a x )X + (a 3 - a 2 )y] 
L*(dy) = a 2 dY + Yd[(a 3 - a x )X + (a 3 - a 2 )Y] . 

It can be checked that L can be extended to the whole sphere. Indeed, it is 
the conformal Killing tensor associated with the spheroconical coordinates 
0. They are the eigenvalues (Ai, A 2 ) of L, since 

det(A7 — L) = A 2 — A[(a 3 - a^X + (a 3 - a 2 )y + (a x + a 2 )] 
+ a 2 (a 3 - ai)X + ai(a 3 - a 2 )Y + a x a 2 
= A 2 (a; 2 + 2/ 2 + ^) 

- A[(a 3 - ai)x 2 + (a 3 - a 2 )y 2 + (ai + a 2 )(x 2 + y 2 + z 2 )} 



a 2 )y 2 + aia 2 (x 2 + y 2 + z 2 ) 



+ a 2 (a 3 — ai)x 2 + ai(a 3 
X 2 (x 2 + y 2 + z 2 ) - A[(a 2 + a 3 )x 2 

+ [a 2 a 3 x 2 + a\a 3 y 2 + aia 2 2 2 ] 
(A - a 2 )(A - a 3 )x 2 + (A - a 1 )(X - a 3 )y 2 + (A - a 1 )(X - a 2 )z 2 



(ai + a 3 )y + (ai + a 2 )z 



x 



+ 



+ 



A — a\ A — a 2 A — a 3 



(A - ai)(A - a 2 )(A - a 3 ) 



Hence, 



(A-Ax)(A-A 2 ) 



x 



+ 



+ 



(A — ai)(A — a 2 )(A — a 3 ) A — a± ' X — a 2 ' X — a 3 

meaning that the eigenvalues of L are the spheroconical cooordinates on 
S* 2 . We have thus deduced the usual separation coordinates of the Neumann 
system simply using the separability condition ([!]) discussed in Section 2. 



5 Integrals of motion for the Neumann sys- 
tem 

The aim of this section is to show that in our setting one can easily de- 
scribe other interesting features of the Neumann system, such as a simple 
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construction of the integrals of motion and the existence of a bi-Hamiltonian 
formulation in a suitable extension of the phase space. The crucial point is 
that the Hamiltonian H of the Neumann system and the complete lift iV of 
the tensor field (§) satisfy the stronger condition 

d(N*dH) = dp! A dH , where p 1 = tr L = X 1 + A 2 , (9) 

implying the separability condition (|l|). This fact has important conse- 
quences, that we are going to show in the general setting considered in Sec- 
tions 2 and 3. 

Let L be a torsionless tensor field of type (1,1), with functionally inde- 
pendent eigenvalues, on an n-dimensional manifold Q, so that we can endow 
the symplectic manifold (T*Q, {•, •}) with the additional Poisson bracket 

{F,G}' = u(NX F ,X G ) , 

using the complete lift iV of L. Let 

det(AJ — L) = X n — { Pl \ n ~ l + p 2 A n - 2 + • • - + p n ) 

be the characteristic polynomial of L, and suppose that H e C°°(T*Q) 
satisfies 

d(N*dH) = d Pl AdH . (10) 

Then one can easily show that ujh — d(N*dH) vanishes on the distribution 
Dh-i so that the separability condition ([!]) is fulfilled and H is separable in 
the coordinates associated with L. We know from Section 3 that there are 
(local) integrals of motion for Xh- If fllOD holds, they can be easily found 
by writing it in the form d(N*dH — p\dH) = 0, and choosing H 2 such that 
dH 2 = N* dH —pidH . Then (as shown in in the case where H is quadratic 
in the momenta) also the 1-form N*dH2 — p2dH is closed, and the process 
can be iterated, so that the integrals of motion for Xh can be determined 
from the recursion relations 

dH 2 = N*dH — pidH 
dH 3 = N*dH 2 - p 2 dH 

(11) 

dH n = N*dH n ^ - p n ^dH . 
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Moreover, we have that 



= N*dH n -p n dH . (12) 

For the Neumann system the global existence of the function H 2 is garan- 
teed by the fact that T*S 2 is simply connected. One finds that, in local 
coordinates, 

H 2 = 2(a lPY 2 Y + a 2 p x 2 X)(X + Y-1)- 2a 3 XY(p Y - Vx f 

1 1 

- -o 2 (oi - a 3 )X - -ai(a 2 - a 3 )Y , 

which coincides with the constant of motion provided by the Lax matrix. 
Notice that the Hamiltonian vector field Xh is not bi-Hamiltonian on T*Q. 
Nevertheless, from (|TJ) we can conclude that 

X H = {;H} = —{.,H n }'. 

Pn 

In the terminology of ||, this means that Xh is a quasi-bi-Hamiltonian vector 
field (see also M, PJ). Next we will show that we can obtain a bi-Hamiltonian 
representation of Xh on the extended phase space M. :=T*Qxl, following 
[18, 0| . The first step is the extension of the bi-Hamiltonian structure from 



T*Q to M.. If c e R is a coordinate in the "additional dimension" and 
F G C°°(T*Q), then the new Poisson brackets are defined as 

{F,c} M = 

{F,c}' M = {F,H}-c{F, Pl } . 

They endow A4 with a bi-Hamiltonian structure. The second step is to notice 
that the recursion relations (|ll|) on the Hi become the usual Lenard relations 
on the functions Hi defined on AA. as 

H = c , Hi = H — cpi for i — 1, . . . ,n. 

Indeed, we have 

{H ,-}m = 
{Hi,-}m — {Ho,-}'m 

{Hu,-}m = {Hn-irjM 
= {H n , -}' M . 
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In the Neumann case, the restriction to c = of the bi-Hamiltonian vec- 
tor field {Hi, -}m — {Hqi'Ym * s ^ ne Neumann vector field. Thus we have 
shown that the Neumann system admits a bi-Hamiltonian formulation in an 
extended phase space, recovering in this way a result of ||. 

6 Concluding remarks 

1. The extension process of a chain of the form (|ll]) into a Lenard chain is 
the opposite of the reduction technique presented in [II], where a suitable 
quotienting produces a chain like flTT| ) from a Lenard chain. 

2. All the results presented here can be extended to the n- dimensional Neu- 
mann system. It would be interesting to compare our approach with the one 
based on the Lax representation. 

3. The condition @ appears also in 0, in the case where {Q,g) is a Rie- 
mannian manifold and H = {<i)ViP j + V(o)- The authors show that (@) 
implies that L is a conformal Killing tensor of g with vanishing torsion. No- 
tice however that in our approach to separability the Riemannian structure 
of Q plays no role and the important objects are bi-Lagrangian foliations on 
bi-Hamiltonian manifolds. 

4. The idea of extending the phase space in order to obtain a bi-Hamiltonian 
formulation of a given (Hamiltonian) vector field goes back, to the best of our 
knowledge, to |Z7|, where it has been used also in the context of separation 
of variables. 
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